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Abstract: We introduce a large class of scalar-tensor models with interactions containing the 
second derivatives of the scalar field but not leading to additional degrees of freedom. These 



<^> models exhibit peculiar features, such as an essential mixing of scalar and tensor kinetic terms, 

which we have named kinetic braiding. This braiding causes the scalar stress tensor to deviate 
from the perfect-fluid form. Cosmology in these models possesses a rich phenomenology, even 
in the limit where the scalar is an exact Goldstone boson. Generically, there are attractor 
solutions where the scalar monitors the behaviour of external matter. Because of the kinetic 
. ^ braiding, the position of the attractor depends both on the form of the Lagrangian and on 

^ the external energy density. The late-time asymptotic of these cosmologies is a de Sitter 

state. The scalar can exhibit phantom behaviour and is able to cross the phantom divide 
with neither ghosts nor gradient instabilities. These features provide a new class of models for 
Dark Energy. As an example, we study in detail a simple one-parameter model. The possible 
observational signatures of this model include a sizeable Early Dark Energy and a specific 
equation of state evolving into the final de-Sitter state from a healthy phantom regime. 
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1. Introduction, Summary and Future Directions 

Infrared modifications of gravity have received considerable attention recently for the poten- 
tially natural explanation of the current cosmic acceleration which they might offer. Much 
advance has been accomplished, with several consistent models now well understood, such 
as (the normal branch of) DGP-type models [1, 2] as well as a number of Lorentz- violating 
setups [3, 4, 5], see also review [6]. An important and intimately related spinoff from these 
investigations has been the study of consistent higher derivative interactions for scalar fields 
[7, 8, 9]. In many modified gravity models the graviton is massive, which implies the presence 
of additional polarizations, the most important of which is a scalar longitudinal mode. This 
polarization can be viewed as the 'Goldstone of a Goldstone' [10], and as a consequence the 
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action for this mode tends to involve its second derivative. This is at the heart of the inher- 
ent difficulty with massive gravity models, since second-derivative interactions easily lead to 
a higher derivative equation of motion and, hence, to ghosts at nonlinear level [10, 11, 12]. 
By exploring the decoupling limit of DGP, in [7] (see also [13]) it was realized that the leading 
cubic interaction for the scalar (j> is of the form (d(jy) 2 \D(p. This interaction leads to an equa- 
tion of motion containing exclusively derivatives of second order, and hence to a ghost-free 
implementation of the Vainshtein effect [7, 13, 11, 14, 15, 16]. More recently, two more inter- 
action terms for a scalar field (quadratic and cubic in ddcft) leading to equations of motion 
of second order have been discussed 1 [8], along with their covariantization [9, 20]. See also 
Refs. [21, 22, 23] where the investigation of a consistent theory of massive gravity in the 
decoupling limit eventually leads to the same additional terms in the Lagrangian. 

Our purpose in this article is to abstract these developments in modified gravity and 
to consider the addition of second-derivative couplings as a generalisation of scalar-tensor 
theories. Specifically, we explore a rather large family of theories where the Lagrangian for 
the scalar field takes the form 



with X = \g^ v d^d v (\> and K,G some generic functions, with K proposed already as k- 
essence [24, 25, 26]. We will show that for any form of the functions K, G the equations of 
motion are of second order. 2 

Hence, this Lagrangian can encode perfectly valid interactions for a single scalar degree 
of freedom. As we will see, a theory like (1.1) has rather unfamiliar properties which we will 
elucidate. We will concentrate on applying such constructions to model Dark Energy. 

The Lagrangian (1.1) can be viewed as reminiscent of a partial resummation of a gradient 
expansion commonly used in fluid dynamics. As is well known, the stress tensor for k- essence 
(the G = case) takes the perfect-fluid form. As we will elaborate on in a future work [27], 
the GO(j) interactions give rise to deviations of a certain type from the perfect-fluid picture. 
Hence, we will refer to the scalar field 4> described by (1.1) as an imperfect scalar. 

Let us also remark on another interesting property of the models of type (1.1). Once 
gravity is introduced, the interaction term GO0 contains a scalar-tensor kinetic coupling, 
schematically G dgd(p, due to the Christoffel symbols present in the covariant D'Alembertian 
operator. A consequence of these couplings is a novel kind of mixing which we have named 
kinetic braiding. It implies that the scalar equation of motion contains the second derivative 
of the metric — and vice-versa — in an essential way: there exists no Einstein frame where the 

1 Similar models were introduced before in another context [17, 18]. Also, recently, such models have been 
extended to vectors and p-forms [19]. 

2 Generic choices of K,G give equations of motion with explicit dependence on the first derivative V M 0. 
Hence, the Lagrangians (1.1) do not have the Galilean symmetry which characterizes the interactions of [7] 
and [8]. 
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kinetic terms are diagonalised 3 . This kinetic braiding manifests itself either at nonlinear level 
or on non-trivial backgrounds. In this sense, this leads to a novel modification of gravity that 
perhaps has not been fully appreciated. Another consequence is the fact that this operator 
may be non-zero and actually important on backgrounds where the second derivatives vanish, 
e.g. the sound speed can be non-vanishing in an exact de-Sitter spacetime. We will give some 
examples of these in the text. 

An obvious domain in which models of type (1.1) can have interesting applications is 
cosmology. As a consequence of kinetic braiding, the Friedmann equation has remarkable 
properties, such as the presence of a term linear in H, similarly to the DGP model. Moreover, 
even in the absence of any direct coupling between the <j> and external matter in the action, 
the evolution of the scalar <j)(t) generically monitors any form of external energy density p e xt 
present. In the shift-symmetric case, where eft is a Goldstone boson, there exist attractor 
solutions on which the shift charge exactly vanishes such that the energy density stored in 
the scalar is a certain model-dependent function of the external energy density, p e ^t- P<f> = 
p^(pext)- Furthermore, for appropriately chosen K, G (whenever |Oj,(0) > 0), the attractor 
tends to de Sitter space in the asymptotic future. This is similar to ghost condensation [3], 
the difference being that on the attractor the ghost condensate behaves like a cosmological 
constant. On the other hand, the imperfect scalar has a nontrivial equation of state dependent 
on the external matter content. 

The existence of this monitoring behaviour opens up a number of applications, the most 
obvious of which is that <f> can play the role of dark energy. Under very general circumstances, 
the on its attractor (j> behaves as a phantom, w* x < — 1. We will show that when evolution is 
started away from the attractor, the imperfect scalar initially redshifts in the usual manner. 
Once the attractor is approached, 4> crosses the phantom divide and ends up dominating and 
accelerating the universe at late times. All of this can be accomplished without generating 
ghosts or gradient instabilities. 

The existence of a physically plausible theory that can cross the phantom divide, i.e. 
violate the Null Energy Condition, opens a Pandora's box full of exotic phenomena such as 
wormholes, bouncing universes, etc. However, currently, wx < — 1 is not observationally 
excluded, see e.g. [28]. 

We analyse in detail a simple one-parameter example that exhibits the features above, 
where the imperfect scalar on the attractor evolves from a subdominant phantom (with wx = 
—2) during matter domination to a dominating late-time de Sitter phase. In addition, the 
energy density of this fluid while off-attractor, has a local maximum at matter-radiation 
equality, allowing us to propose this model as a simple example of a form of early dark 
energy. 

This article is organized as follows: We introduce the class of models (1.1) in section §2, 

3 An ordinary scalar with nonminimal coupling to curvature, e.g. (j> 2 R, also leads to kinetic mixing. However, 
in this case the mixing is non-essential in the sense that it can be removed by a field redefinition (in the form 
of a conformal transformation). Note that kinetic braiding is present in the higher-order terms of the covariant 
Galileon [9, 20]. 
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discussing their equations of motion and coupling to gravity along with the effective metric for 
the propagation of perturbations. We specialise our results to the cosmological background 
in section §3, and discuss the stability conditions for background solutions: positivity of the 
sound speed and of the kinetic term of perturbations. In section §4, we specialise to the 
cosmology of the shift-symmetric case. We discuss the dynamical attractors in the presence 
of matter and in the late universe — once matter has diluted away. In section §5, we discuss 
a simple one-parameter model to demonstrate explicitly the sort of phenomenology one can 
obtain within the new framework proposed. We invite those readers more interested in dark- 
energy phenomenology to start there. Finally, in Appedix A, we present a derivation of the 
quadratic action for cosmological perturbations for theories described by (1.1). 

This work has a certain overlap with Refs [29, 30, 31, 32, 33] where the cosmology of the 
decoupling limit of DGP and some of its simple modifications are considered. These references 
concentrate on models with G = Xf(4>) and non-minimally coupled to the Ricci scalar, a la 
Brans-Dicke gravity. In this work, we extend the above class of theories to a general function 
G(4>, X) while restricting to a minimal coupling to gravity. As such we make it clear that the 
new interesting phenomenology is a result of kinetic braiding. In [34, 35] the authors take the 
effective field theory approach to discuss the ability of generalisation of the decoupling limit 
of DGP to stably violate the Null Energy Condition and propose a new mechanism to replace 
inflation as the initial stage of the universe. In addition, in Refs [36, 37], the cosmology of 
higher-order terms of the general covariant galileon is analysed. 

We believe that the models described by Eq. (1.1) provide many avenues for exploration. 
In particular, the reinterpretation of these models as imperfect fluids yields many interesting 
insights, as we will show in our upcoming work, Ref. [27]. It is well known that the classical 
results of single- field inflation can be obtained entirely in terms of perfect-fluid variables [38] . 
The implementation of inflation using imperfect scalars provides another rich arena for model- 
building, with potentially non-trivial effects on cosmological perturbations. We will return to 
a thorough investigation of cosmological perturbations in the future. 

2. General Properties of Kinetic Gravity Braiding Theories 

In this paper we consider a class of scalar field theories minimally coupled to gravity, which 
is described by the action 4 



and K (<f>, X) and G (</>, X) are arbitrary functions of the scalar field and its standard ki- 
netic term X. Here and in most of the paper, we use reduced Planck units where Mp\ = 

4 We use the metric signature convention (H — ). 





where we denote 




9' 
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(87tGn) = 1. In further discussion we will also use the corresponding Lagrangian for the 
scalar 

£(4>,X,B) = K((f>,X) + G((f>,X) xB, 

which we will consider as a function of three variables: (/>, X and B; whenever we differentiate 
with respect to <j) or X, we are keeping B constant. Our class of models extends those studied 
up to now by generalising the interaction term X of DGP decoupling limit and the galileon 
models to a function G(<f),X), significantly enriching the phenomenology. 

One can naively motivate our interest in non-canonical field theories containing the new 
G x B term by a similar argument to that for k- essence. In k- essence, we treat the canonical 
kinetic term X as a first term in an effective-field-theory expansion, since no symmetries exist 
that would prevent additional contributions, i.e. 

K((f),X) ~ X (l + Cl (<j))X + c 2 (<f))X 2 + ...) . 

Then after resummation of these corrections one obtains a non-canonical Lagrangian with 
a derivatively self-coupled scalar field, e.g. as in the case of the Dirac-Born-Infeld theories. 
However, one could equally well start from the canonical kinetic term X and integrate it by 
parts with no effect on the observed dynamics. In that case, it would not be unreasonable to 
expect there to exist higher-order corrections similar to those discussed above, 

G ((f), X) Dcf> ~ -</n<t> (1 + ci ((f)) X + c 2 ((f)) X 2 + ...) . 

Note that both these expansions are in X. Below we will show that it is not possible to recast 
such corrections to —(f>\Z\(f> as a modification of the function K ((f), X) and surface terms in the 
action, and therefore these theories are outside of the k- essence scheme. We stress that, in the 
Einstein-Hilbert action, the second derivative enters the Lagrangian only linearly and, despite 
the presence of these higher derivatives in the Lagrangian, the equations of motion remain of 
the second order. Therefore, similarly to the Einstein's general relativity, there are no new 
(ghosty) degrees of freedom which would appear as the result of the Ostrogradsky theorem 5 . 
Thus, at least for some functions K ((f), X) and G ((f), X) and some general backgrounds, there 
are no ghosts and pathologies associated with them — the theories containing G x B are on 
the same footing as k- essence or DBI theories. 

The form of the action in Eq. (2.1) can be further generalised by introducing a non- 
minimal coupling of the scalar field to gravity, or a modification of the f(R) type. An 
appropriate change of variables will transform the action to the Einstein frame, where the 
usual Einstein-Hilbert term is restored. In the case of a theory containing only gravity and the 
scalar, this is equivalent to choosing different functions K and G. However, if external matter 
is present, this transformation introduces an additional coupling between the matter and the 

5 There is another important similarity with General Relativity: a rigorous formulation of the variational 
principle requires a boundary term similar to the Gibbons-Hawking- York term in General Relativity. In this 
work we do not consider this boundary term assuming that it exists. The term appropriate for the decoupling 
limit of DGP was found in [39] . 
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scalar field 6 . Such theories can actually evade Solar-System tests of general relativity because 
of the presence of the Vainshtein mechanism [14, 47]. However, in this work we will concentrate 
on the minimally coupled models. This minimal coupling to gravity is the crucial difference of 
our work from the recent studies of the Galileon scalar-tensor models [29, 30, 31, 32, 33, 48]. 
In particular, we have to stress that the most unexpected and interesting features of these 
theories are a result not only of the possible non-minimal coupling of the scalar to gravity 
but of the presence of G x B in the action. 

In our discussion we will frequently use a Lagrangian equivalent to C: integrating by 
parts the scalar- field contribution to the action Eq. (2.1), we obtain 



V = K 



V x 6) V 



G = K - 2XGs - G X VVV A X , (2.2) 



where the subscripts <j> and X denote partial differentiation with respect to these independent 
variables. It is clear that it is the dependence of G on the field's gradient, X, that prevents 
this theory from being recast as a k- essence model, whereas the derivative of G with respect 
to can be recast as an additional k- essence term. We will therefore only consider such 
non-trivial models where Gx ^ 0. The above form of the Lagrangian also implies that no 
quantities will ever be dependent on the function G (<f>, X) itself, but only on its derivatives. 

2.1 Equation of Motion for the Scalar Field 

The equation of motion for the scalar field can be obtained in the usual way by varying the 
action 

= ^ - V„ ((Cx - G+) W - VG) , 



v-g d( t 

where V$ = K§ — [(V A t/>) Va] G^. In particular, this presentation of the result makes explicit 
the existence of a Noether current for Lagrangians invariant under constant shifts of the field, 

4> — > 4> + c, 

J M = (C x - 2G4,) - GxV^X . (2.3) 
In terms of this current, the equation of motion takes the elegant form 

V M J^ = V* . (2.4) 

The fully expanded equation of motion is somewhat unwieldy, but its structure is important 
for future discussion: 

L^V M V^ + (VaV^) Q a ^ u (V M V^) + - G x R^V^V u <f> = (2.5) 

6 Following the same logic one could also consider additional derivative couplings of the scalar field to the 
Einstein tensor. In that case, equations of motion would still be of the second order, see e.g. recent works 
[40, 41, 42, 43, 44], contrary to [45, 46]. 
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where R^ v is the Ricci tensor and £, L^ u and Qap^v are constructed out of the metric g^ v , 
field <f> and its first derivatives V^^> only. The term 8^ which is zero order in the second 
derivatives is the field derivative of 

£ = 2X(K X -G (f) )-K, (2.6) 
the terms linear in the second derivatives are contracted with 

L„ u = (K x - 2G + 2XG X(t >) + (K XX - IG^x) V^V„tf , 
and the terms quadratic in the second derivatives are contracted with 

QaPnu = 1 (g<xfl H ^ + H aPgi^ _ I [g^H ua + g^H^ + g>* a H v P + g ua H^ , (2.7) 
where 

Hp, = G x g^ + GxxV^V^ . (2.8) 

Most importantly, the equation of motion contains at most second-order derivatives, and 
therefore the theories with kinetic braiding do not contain any hidden degrees of freedom. 
Moreover, the structure of the tensor Q a ^^ u implies that <ft appears only linearly in the 
equation of motion and therefore the equation of motion is solvable with respect to (f>, for 
any time coordinate 7 . This is sufficient to ensure that this equation is of the normal type 
required by the Cauchy-Kowalewski theorem and that the Cauchy problem has a unique local 
solution, at least for analytic functions. 

The presence of a third-order derivative might have been expected owing to the d'Alamber- 
tian term in the Lagrangian; however, these terms arise in the combination (VD — DV) <f> and 
can be commuted away, leaving behind a term coupling the scalar to the Ricci tensor, R^. 
This is the key feature of theories with kinetic braiding. Namely a seemingly minimally cou- 
pled theory has an equation of motion depending explicitly on the curvature so that the scalar 
equation of motion includes the second derivatives of both the scalar field and the metric. 

It is interesting to note that Eq. (2.5) can be considered as a generalisation of the Ampere- 
Monge equation for four-dimensional manifolds with Lorentzian signature. 

2.2 Energy-Momentum Tensor 

The energy-momentum tensor ("EMT") for the scalar field is most easily derived in the 
standard way using the Lagrangian presented in Eq. (2.2) 

2 SS 

Tfv, = —= = CxV^V^ - g^V - - V U GV^ . (2.9) 

It is key to observe that the EMT also contains second derivatives of <f> appearing in V, £x 
and V^G. Therefore both the Einstein's equations and the equation of motion for the scalar 

7 Lorentz symmetry implies that this statement is true for the second derivatives <j) xx with respect to any 
coordinate x 
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field, Eq. (2.5), contain second derivatives of both the metric g M „ and the scalar field (j>, so 
that the system is not diagonal in second derivatives. This kinetic mixing is essential and 
cannot be undone through a conformal transformation. We shall refer to such a mixing as 
kinetic braiding and show that it leads to a number of surprising properties for this system. 

In particular, for general configurations with timelike field derivatives the structure of the 
EMT is not that of a perfect fluid, contrary to minimally coupled k- essence theories, including 
canonical scalar fields. Moreover, for some of such configurations, the EMT does not even 
have a timelike eigenvector. The imperfect nature of the fluid is directly related to kinetic 
braiding. In the paper [27], we discuss the corresponding imperfect-fluid picture in detail. For 
cosmological solutions, which are exactly homogeneous and isotropic, the EMT is forced to be 
of perfect-fluid type. However, cosmological solutions including small perturbations possess 
an EMT with small terms deviating from the perfect fluid. This changes the standard picture 
of cosmological perturbations where the intuition was gained from the perfect fluid case. 

Using the Einstein equations including the EMT contribution T^f of fields different from 

Rnu — -zg t j,uR = T l j iV + T^ V , (2-10) 

we can eliminate the second derivatives of the metric in the scalar field equation of motion 
Eq. (2.5). In this way we obtain 

irv^vt + (v Q v^) Q a ^ v (v M v v 0) + z- (2.11) 

-G x - 1<TtJ\ V^Vut = , 

where the new term which does not contain second derivatives is 

Z = S 4> - 2XG X (X (K x - 4G^) + K) , 

while the terms linear in the second derivatives are now contracted with 

Lya, = V - 2XG\ (Xg^ - 2V^V„0) . (2.12) 

The tensor Q a ^^ v remains the same. Note that the new contributions modifying £ and 
are suppressed by M™. In Eq. (2.11) we have obtained a form of equations of motion which 
contains the second derivatives only of the scalar field <fi but not of the metric 8 . This is im- 
portant for the investigation of causality and stability in this class of models, see section §2.3. 

Let us finally emphasise the fact that the scalar inevitably becomes coupled to external 
matter in the non-linear fashion resulting from the last term of Eq. (2.11). This occurs 
in a quite peculiar manner: unless an additional explicit direct scalar/matter interaction is 
introduced, this coupling is one way and external matter does not become coupled to the 
scalar, only feeling it through its gravitational effects. This coupling will play a central 
role in the following sections, is a consequence of the kinetic braiding. A proper analysis of 

8 Here we assume that, as usual, other matter fields do not contain second derivatives in their EMT. 
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the constraints on the model arising as a result of this coupling will be deferred to future 
investigation, but for some headway in this direction see Ref. [33]. Let us only remark here 
that the self-interactions of the scalar can give rise to a Vainshtein mechanism that helps to 
evade the observational constraints. 

2.3 Effective Metric 

Let us perturb the system comprising the scalar equation of motion and Einstein equations 
around a given solution <p{x^) and g a s (x^). Without loss of generality, we will neglect the 
perturbations in the external fields, so that 5T®„ = 0. Now let us find the characteristic 
surfaces S (wave-fronts) for the linearized system. This is crucial for the investigation of 
causality and stability with respect to the high-frequency perturbations, see e.g. [49]. 

The linearized equation of motion for the scalar field Eq. (2.11) does not contain sec- 
ond derivatives of the metric. Therefore the matrix of the principal part of the differential 
operator for this system is lower triangular in {5(f), 5g^ u ) coordinates. This implies that the 
characteristic equation for the system (see e.g. [49, p. 580]) becomes just a product of a 
purely gravitational and a purely scalar part and the second derivatives of the scalar field 
still present in the linearised Einstein equations do not play any role for characteristics. The 
causal structure for the propagation of the scalar perturbations can therefore be directly read 
from just equation (2.11). Note that the same procedure for the original equation of motion 
(2.5) would give an incorrect result. 

From the standard eikonal ansatz 5(f) = A {x) exp {iuiS {x)) with a slowly varying ampli- 
tude A {x), taking the formal limit co — > oo, we obtain 

g" v d^Sd v S = 0, (2.13) 

where the effective contravariant metric for propagation of perturbations is 

guv = l»v + 2Q a ^ v V a V p( f) . (2.14) 

Using formulae (2.7) and (2.12), we can expand this expression as 

9w = ^g„u + - (G x V v <f>) - Vu (GxV^) , (2.15) 

where we have introduced the notation 

Q = C XX + AXG 2 X , (2.16) 

and 

fi = C x - 2G^ + V A (G*VV) - 2X 2 G\ . (2.17) 

The velocity of the wavefront S and the corresponding causal structure of the acoustic space- 
time should be inferred from the inverse or cov&ricuit metric Q^j} for which Q ^\ Q^ v — 5^, 
for details see e.g. Appendix A of Ref. [50]. The effective metric for the decoupling limit of 
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the DGP model (where G(4>,X) oc X) was first obtained in Ref. [51], see also e.g. [34]. The 
result (2.15) is more general, because it is derived for a generic function G (cp,X). 

The main obstacle for the analysis of this metric for general backgrounds is that Q 
contains not only second derivatives in f2 and but also an additional tensorial structure 
which is a Lie derivative of the metric £yg iiv = V M (Gx^v4>) + V v (GxV M </>) with respect 
to = Gx V^</>. This structure prevents V u<j> from being an eigenvector of the metric for a 
general background as was the case for k- essence. 

In this paper, we will restrict our analysis to the cosmological case. In particular, in 
section §3.2, we calculate the sound speed and present the condition for the absence of ghosts 
using the metric (2.15). The same results are confirmed by a direct calculation of the action 
for cosmological perturbations, presented in Appendix A. The advantage of this effective- 
metric formalism is that it allows for the study of causality and stability for high-frequency 
perturbations for any given background without calculating the action for perturbations, 
which may be a very complicated task for general, less symmetric backgrounds. 

3. Cosmology 

3.1 Background Evolution 

In this section, we will discuss cosmology in braided models (1.1). We will first consider the 
homogeneous and isotropic background solutions 4>{t) for a generic model and, in Section 4, 
we will specialize to the shift-symmetric case, where the integration of the equations of motion 
is almost straightforward. 

For simplicity, we will restrict our discussion to the spatially flat Friedman-Robertson- 
Walker ( "FRW" ) metric 

ds 2 = dt 2 -a 2 {t)dx 2 . (3.1) 

On this background, we have X = \(f> 2 , where the dot represents the time derivative. In the 
current Eq. (2.3) only the "charge density" 

J = J = (k x - 2G <i> + 3H<j>G x ) <j> , (3.2) 

survives, so that the equation of motion (2.4) takes the form 

J + 3HJ = Vtj, , (3.3) 

where, as usual, H = a/a denotes the Hubble parameter. Note that this equation of motion 
not only contains <j) but also a. As a result of the symmetry of the cosmological setup, the 
energy-momentum tensor takes the perfect-fluid form and the only non- vanishing stress-tensor 
components, Eqs (2.9), reduce to the energy density 

e = T° = 4>J + 2XG<j> - K , (3.4) 

and the pressure 

V = -^T\ = K-2XG lj} -2XGx4>, (3.5) 
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which turns out to be given by the Lagrangian (2.2). The two unusual properties arising as 
a result of kinetic braiding are that the pressure contains 4> while the energy density depends 
on, H, the Hubble parameter 9 . 

It is convenient to introduce the variable "measuring" the strength of the kinetic braiding 

k = XG X ; (3.6) 

its physical meaning will be further elucidated in the paper [27]. 

The part of energy density (3.4) which is not braided and which does not depend on the 
Hubble parameter looks similar to k- essence: 

£ = 2X(K X -G^)-K, (3.7) 

and its field derivative appears already in the equation of motion (2.5), so that 

e(0, <j>, H) = £ + 6K<pH . (3.8) 

The Friedmann equation becomes not the usual complete square but a general quadratic 
expression for H, 

H 2 = 2k<PH +\{£ + Pcxt) , (3.9) 

where we have introduced the energy density p cx t from any additional sector external to 
(f). This could be dust, radiation, a bare cosmological constant term or any other standard 
cosmological fluid. 

It is useful to compare this equation with the Friedmann equation arising in the DGP 
scenario [53, 54]. In particular, we see that {n(f))~ l plays the role of the crossover scale rn 
beyond which gravity starts probing the extra dimension in DGP, while £ corresponds to the 
energy density localized on the brane. 

We can solve for H in Eq. (3.9) to obtain 

H = K<j> + aJ (4) 2 + ~(S + pext) , (3.10) 

with a = ±1 denoting a choice of branch. We will see in section §4.3 that only one choice of 
a (which depends on the forms of K, G) is compatible with stable fluctuations around the 
solution. 

The presence of the term linear in H in (3.9) is responsible for a mismatch between the 
two branches of solutions (labelled by a) and expansion or contraction of the cosmologies. 
Only a simultaneous transformation (a = +1,4') — > (cr = gives the time-reserved 

solution. 

For future reference, it is useful to also write down the acceleration equation, 

H = --(e + /9 e xt + V +Pcxt) = K(j>- -J(f> - ^(Pcxt +Pcxt) , ( 3 -ll) 

9 Note that this situation is different from the so-called Inhomogeneous Equation of State where pressure is 
postulated to be a function of both the energy density and H [52]. 
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where p ex t is the pressure of the external matter; we will also sometimes use its equation-of- 
state parameter, w ext = Pext/Pcxt- 

Finally we present the expanded form of the equation of motion (2.11) in a cosmological 
context 

Dip + 3J(H - 4) + e^ = 3K(p cxt + p ext ) , (3.12) 

where 

D = Ex + 64>Hkx + 6k 2 . (3.13) 

In the absence of kinetic braiding, k = 0, this reduces to the ^-essence equation of motion, 
c.f. Ref. [55, Eq. (11.13)]. One notices that the external matter content explicitly filtered 
down into Eq. (3.12) despite the absence of any direct coupling between (ft and matter in the 
Lagrangian. This is traced back to the explicit dependence of J on H and is a consequence 
of the braiding between the scalar and gravity in these models. The normalization factor D 
in Eq. (3.13) will appear repeatedly below and it plays a prominent role since it determines 
whether the fluctuations of (ft around the background are ghosty (D < 0) or not (D > 0), see 
section §3.2. 

Notice the presence of a term quadratic in k. Once the Mp\ factors are restored, we see 
that this term is suppressed by one more power of Mpj, reflecting that it is a type of backre- 
action effect implied by the scalar-gravity braiding. This term arises from the dependence of 
the equation of motion for (ft, Eq. (2.5), on H and the simultaneous dependence of H on (ft as 
implied by the Friedmann equation. 

In the absence of external matter, close to the point where J = 0, (ft will be small provided 
e is a slowly- varying function of (ft, through Eq. (3.12). Then Eq. (3.11) implies that we are 
dealing with a generalised slow-roll, quasi-de-Sitter state. 

3.2 Stability 

Now let us apply the machinery of the effective metric to cosmological solutions. Plugging in 
the cosmological background into the general expression for the metric (2.15) we obtain for 
the contravariant metric 

= 5%5%D - ^b v ia ~ 2 (n - jH^j . (3.14) 

From this expression, it follows that ghosts are absent if D > 0. Indeed, it is Q^ v which 
is contracting the derivatives in the kinetic term for perturbations 10 , see (A. 8). Here it is 
important that the sign in front of be chosen correctly — so that this condition continuously 
transforms to that of k- essence. 

10 In the high wave-number limit, the gauge-invariant canonical variable v is proportional to the gauge in- 
variant 5(j>. Also note that the full calculation of the quadratic action gives the expression for the normalisation 
of the cosmological perturbations, Eq. (A. 9). This is modified compared to the standard case and therefore is 
likely impact predictions of inflation. 
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This contravariant metric is diagonal. Therefore it is easy to find its inverse (or covariant) 
metric along with the corresponding acoustic line element 

dS 2 = g^dxW =(n- i?-Hj (c 2 s dt 2 - a 2 (t) dx 2 ) , (3.15) 



where the sound speed is given by 

2 Q.(j> — 4kH 



(3.16) 



Thus for the stability with respect to high-frequency perturbations we should require that 11 

4k 

D > and - -^H > 0. (3.17) 

Only those models for which during the whole evolution history these two conditions are 
satisfied can be considered as plausible from a physical point of view. Further we calculate 
fi for the cosmological background and obtain 



dK J + 2 ( k + Hk — K 
n-^H= ^ . (3.18) 



while the definition (3.13) can be rewritten as D<j) = — QHk + 6k 2 where the partial 
derivative of the energy density (which is considered as a function of three independent 
variables, see Eq. (3.8)) is taken keeping H = const, as if the Hubble parameter were externally 
fixed. 

Using these expressions, the formula for the sound speed (3.16) can be written as 



J + 2k + 2k I H - Kq 
c 2 = ^ . (3.19) 



e ; — 6k 



(H - K<ty 



This expression depends on 4> which can be eliminated using the equation of motion (3.12). 
It is important to note that 

-p ■ + 4k + 2k (AH - k 



Cs = — —A ry^- + ~ , (3-20) 

e ; - 6k [ H - k6\ e 



even in the shift-symmetric case, contrary to k- essence [38]. Here the derivative of the total 
pressure is taken keeping eft = const. The sound speed in these models is arbitrary and can 
even be larger than the speed of light. In particular, this is exactly the case in the model 



11 Note that the tensor modes are not affected in our model, as can be seen in the full calculation of the 
quadratic action Eq. (A. 8), and therefore do not provide any new conditions for stability. 
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discussed in the section §5, see e.g. Fig. 4. In this respect the situation is not very different 
from k-essence 12 . From the expression for the acoustic interval Eq. (3.15), it follows that 
there are no causal pathologies (Closed Causal Curves) even in the presence of super luminal 
propagation. We follow the discussion of this issue presented in Ref. [50] (see also the most 
recent paper [57] and older works [58, 59, 60, 61]) and claim that this superluminality does not 
imply any inconsistencies. A different issue is a possible UV completion of the theories with 
superluminal propagation. Currently, there are good arguments [62] that this UV completion 
cannot be realized within a Lorentz-invariant and renormalizable local QFT or perturbative 
string theory framework. However, on the level of an effective-field-theory description there 
are no inconsistencies. This is also accepted by some of the authors of [62], see [34, 56]. This 
possibility to realise superluminal propagation of perturbations on nontrivial backgrounds 
can lead to interesting cosmological scenarios in the context of inflation [63] and some exotic 
alternatives, e.g. [34, 64, 65, 66]. Moreover, this superluminality raises important questions 
regarding black hole physics [67, 68, 69, 70]. 

4. Cosmology of Shift-Symmetric Gravity-Braided Models 

In this section, we will restrict our attention to models that realise an exact shift symmetry 



in terms of an appropriate choice of the scalar-field variable (p. In practice, this implies that 
K, G must be (^—independent, that is the scalar Lagrangian takes the form 13 



Imposing such a symmetry, i.e. requiring that the field only be derivatively coupled, means 
that (j> can be interpreted as a Goldstone boson of some broken symmetry. This naturally 
prevents such a field from acquiring a mass. As we will see, the dynamics of such kinetically 
braided Goldstone bosons make them a very compelling model for dark-energy dynamics. In 
addition, the scalar would have to be derivatively coupled to matter, if at all, which would 
make evading fifth- force constraints relatively easy 14 . 

12 However, contrary to the case of k-essence and the conclusions of Ref. [56], in the model presented in 
section §5 there is a period during which the null-energy condition is violated on an isotropic background while 
the speed of sound is sublumninal, see Figs 4 or 5. 

13 Shift-invariance allows for somewhat more freedom in choosing Lagrangians. Obviously, a term in G linear 
in <f> is compatible with this symmetry. However, it can be absorbed in the form of K by integration by parts. 
Similarly, a change of variable <j> — > <j>(4>) w iH trivially introduce </>-dependence but will still realize the same 
symmetry, even though not as a shift in (p. 

14 Generic derivative couplings to matter can still give rise to sizeable effects, but typically the most prob- 
lematic couplings can be suppressed by additional symmetries [3]. In fact, in this sense the situation in the 
present model quite similar to ghost condensation [3] and Hofava gravity [71, 72, 73] 



(p —?■ (f) + const 



C = K(X) + G(X) U<j>. 
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In the shift-symmetric case, the shift-current (2.3) is conserved, and the scalar equation 
of motion is precisely this statement, V M J M = 0. In a homogeneous and isotropic background, 
this reduces to simply 

j + 3HJ = 0, (4.2) 

that is, the shift charge in a comoving volume is constant. Furthermore, the equation of 
motion Eq. (4.2) can be trivially integrated, 

• ■ const 

J = <P(K X + tyHGx) = jjTjT • (4-3) 

Hence, the expansion of the universe drives J to zero, and the locus of J = in configuration 
space represents a future attractor for expanding FRW space-times 15 . As we will see, even 
on the attractor one can have rather interesting behaviour, basically because the value of (p 
on the attractor depends on the amount of external matter density at every moment of time 
(because of the appearance of H in J). 

It is convenient and illuminating to split the energy density of the scalar into the contri- 
bution from the attractor and that of the departure from it, 

e = e* + e J , 
with e* = e|j=o 

and e J = e-e* . (4.4) 

Both 'components' behave like fluids with a particular equation of state on the background. 
By definition, the J (or 'off-attractor') component dilutes away with expansion, so it behaves 
like a rather standard fluid. As we will see, the attractor component can be quite exotic, 
as it can easily exhibit phantom behaviour without generating instabilities. In fact, we will 
find that requiring that the energy density in the attractor be positive and that the field 
fluctuation be healthy (not a ghost) implies that generically the attractor displays phantom 
behaviour. Hence, the composition of the J- and attractor components (that is, the situation 
with a generic initial condition) is such that the energy density in <fi first dilutes away and 
then grows. That is, we will find that <fi can cross the phantom divide without leading to 
instabilities. 

At this point, we note that since the action contains second derivatives, therefore possi- 
bility of smoothly crossing of the wx = — 1 barrier does not contradict the statement proved 
in [55] and rederived in different ways later in [75, 76, 77, 78, 79, 80, 50]. Therefore kinetic 
braiding provides a working example of the so-called Quintom scenario of Ref. [81], see also 
reviews [82, 83, 84]. Thus in this respect, kinetic braiding exhibits similar phenomenology to 
models with explicit nonminimal coupling to gravity, which also allow one to have a classi- 
cally stable crossing of the phantom divide in scalar-tensor theories [85, 86, 87, 88]. Another 

15 The appearance of de-Sitter attractors in shift-symmetric theories with higher derivatives was also used 
in the so-called the B-Innation [74]. This behaviour was first noticed by Alexei Anisimov, one of the authors 
of the B-Infiation paper, who regrettably passed away after this preprint was submitted to the arXiv, RIP. 
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theory comprising a single degree of freedom which is able to penetrate the phantom divide 
without classical instabilities is the so-called A 93- fluid ("Angel Dust") [89]. However, this is a 
very non-standard theory where issues related to the strong-coupling scale and quantisation 
remain to be addressed. In the case of kinetic braiding, the quantisation of perturbations is 
standard and the crossing occurs in a regime free of negative energy perturbations and which 
is completely under control. Models which exhibit healthy phantom behaviour within the 
Galileon framework were realised in Refs [34, 35]. Moreover, in Ref. [90] a healthy violation 
of the Null Energy Condition was achieved by including (d(p) 2 on the level of perturbations 
within the effective-field-theory framework. 

A central issue in the significance of the split (4.4) is if and when the off-attractor compo- 
nent dilutes away sufficiently fast so that the attractor regime is actually reached. We discuss 
this among other points in §4.2. But first, in the next subsection, we shall analyse the generic 
properties of attractor solutions. 

4.1 The Phantom Attractor 

First of all, since the Friedmann equation fixes H as a function of X = (fi 2 /2 and the external 
matter energy density p ex t, one can view J = (j)(Kx + 3<j)HGx) as a function of the two 
variables X and p cxt , 

J = J(X,p ext ) , 

whose form depends only on the choice of K, G. The attractor solutions are the roots of 

J(X„p cxt ) = (4.5) 

We will furnish quantities evaluated on the attractor with an asterisk ('*') sub- /superscript 16 . 
In particular, for any form of external matter, the attractor behaviour of the scalar is such 
that X becomes a certain local function of p e xtj 

X* = X4 Pcxt ), (4.6) 

the form of which is solely determined by the dynamics, i.e. by the form of K, G. Notice that 
this holds irrespective of the time dependence of p ext (i-e. its equation of state), so implicitly 
Eq. (4.6) also represents the time evolution of X(t) on the attractor, X(t) = X*(p ext (t)). 

From (3.4) it follows that on the attractor the energy density stored in the scalar, e*, is 
also only a function of p ex t, 

e* = e*(p ex t) = -K(X.( PfSKk )) . (4.7) 

Hence, the attractor has the remarkable property of responding to the external energy density 
in a way determined by the form of K and G. 

16 In principle, one should introduce a label to distinguish among the several possible attractors each of 
which has its own non-overlapping basin of attraction. In addition, the form of J(X, p ox t) also depends on the 
branch for the Friedman equation (er), so the attractors in general also depend on a. For the sake of clarity, 
we will spare the reader these labels with no risk of introducing ambiguities. 
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We should emphasize that since e* is a function of p ex t only, this effectively allows one 
to modify the Friedman equation (viewed as the relation between H and p ey ±) in almost any 
way in these models. K and G can be picked in such a way so as to replicate any evolution 
history, once the shift current dilutes and the attractor is reached. The attractor present 
in our kinetically braided model can therefore be considered as a concrete realisation of the 
phenomenological "Cardassian" scenario proposed in Ref. [91]. 

Of course, one must make sure that the perturbations are healthy, which will not be the 
case for a generic choice. It should be noted that these modifications are the result of the 
fact that scalar-gravity braiding modifies the Hamiltonian constraint in gravity. 

The behaviour of the attractor can be easily understood in terms of the function e*(p ex t). 
For example, if e*(p e xt) is a decreasing function (such as the one in the model of section §5), 
then as matter dilutes the attractor energy density grows, implying a phantom equation of 
state. In the presence of usual external matter (io ex t > — 1)> e * will eventually dominate the 
expansion of the universe and, as long as e* (0) > 0, the final state will asymptote to de Sitter. 
We will see shortly that it is precisely this kind of setup, with phantom behaviour for the 
attractor, which is free of instabilities. 

At this stage, the form of the function e*(p ex t) is inevitably implicit — below we will 
elaborate on how to extract it in general once K and G are given. However, for most purposes 
its explicit form is not really necessary. 

For instance, the effective equation of state of the attractor can be written as 



Since the normalisation factor D* > 0, so that the <fi fluctuation not be a ghost we find that 
(assuming u> ext > — 1): 

• the attractor is phantom (w x < — 1) whenever its energy density e* = —K* is positive, 
and vice-versa; 

• as the external energy density dilutes away, the attractor approaches de Sitter with 
w* x — > — 1 from below — the phantom region, provided e* > 0; 

• the equation of state of the attractor contribution is tied to the equation of state of the 
external matter and therefore will change, for example, at matter-radiation equality. 

It is quite simple to imagine that it is the behaviour of the attractor that is relevant for 
the observations of dark energy. Therefore, provided that the off-attractor contribution has 
diluted away early enough during the history of the universe, it is a generic prediction of the 





The last equality follows from the chain rule and reading off 
equation of motion on the attractor, 




t-) 7 Kcxt 
<p* — At* 



(4.9) 



models of Eq. (4.1) that dark energy will be seen to violate the NEC and approach de Sitter 
from below — the phantom region. 

Equations (4.7) and (4.8) encode the most important feature of the behaviour of the 
attractor, namely that <f> responds to the presence of any other form of matter, and that a 
stable quasi-de Sitter stage is reached only if the attractor is phantom. So far, we only used 
the positivity of the energy density (— -K*) and -D*. An additional constraint required for 
stability is discussed in §4.3. 

Let us close this subsection with some gymnastics that can ease the task of finding the 
attractor, in particular the function e* (pext)j f° r an arbitrary K and G. The technical difficulty 
is that the explicit expression of J(X, p ex t) found from Eqs (3.2) and (3.10) is rather involved 
and may not allow to solve easily for its roots. However, on the attractor, the equations 
simplify considerably — for example, the Friedman equation becomes 

3H 2 = -K*+p cxt . (4.10) 

From the equation J = itself we see that we can also write = — K x / (3<f)*G x ) . Thus, it 
is easy to see that solving J = is equivalent to solving 17 



+ = p cxt and (4.11) 



(K*-X*K* X \ ■ , A , n . 

sgn ^ — j = o-sgn0* . (4.12) 

Clearly, this form of the equations for the attractor is considerably simpler. In particular, 
Eq. (4.11) directly gives p e xt as a function of X*, which is just the inverse of the function in 
which we are interested, X*(p cx t). Also, it is straightforward to obtain the relation between 
the attractor energy density e* and p e xt m parametric form, 

in terms of the 'parameter' X. This curve in the e*-p ex t plane encodes how the attractor 
responds to the external matter p ex t for given K, G (the a = ± branches corresponding to 
consecutive segments of the curve). 



Finally, let us mention that the roots of 

K 2 



§X*G 2 X ^ 



+ K* = (4.13) 



represent the pure-de-Sitter attractors (when K* < 0), corresponding to the asymptotic state 
of expanding solutions, once all matter and shift-current density, J, dilute away. These are 

17 Eq. (4.11) needs to be supplemented with (4.12) because the l.h.s. of Eq (4.11) contains less information 
than J. Specifically, Eq (4.11) does not involve H and hence is independent of the choice of branch of the 
Friedman equation, a. Eq. (4.12) provides that information. 
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kinetic condensates similar to ghost condensation [3] in that the non-trivial condensate 
mimics a cosmological constant. However, let us emphasize that the present situation is 
qualitatively different from ghost condensation for two reasons: First, in ghost condensation 
(the G = case) the attractor can only behave like a pure cosmological constant. Instead, 
with G ^ it has a nontrivial (phantom) equation of state in the presence of external matter. 
Second, the fluctuations around the attractor behave very differently in ghost condensation 
and the model with kinetic braiding (4.1). Even on the pure de Sitter condensate the sound 
speed squared of the scalar fluctuation can be positive, Cg > 0. This is a consequence of the 
presence of the of the term and the resulting imperfection of the effective fluid. As 

we will see in §4.2, around the de Sitter condensate, the equation of state corresponding to 
a homogeneous perturbation (the off-attractor component due to J ^ 0) vanishes. Hence, 
the off-attractor component behaves like a pressureless fluid. For a perfect fluid this would 
imply that = 0, as indeed happens in ghost condensation (which is why one needs to resort 
to higher space-derivative terms in Ref. [3]). In our case, despite having w J = for the 
off-attractor component, one finds Cg > (as long as H ^ 0). This can also be seen from the 
EFT point of view by including operators like GO</> in the Lagrangian for perturbations [90]. 

4.2 Approach to the Attractor 

The background dynamics away from the attractor are given by Eq. (4.3). The quantities 
that evolve with time in a simple way are J ~ a~ 3 and p e xt ~ a _3 ( 1 + u, ext)_ j n principle, then, 
to explicitly work out the time evolution of the system we need to express everything in terms 
of J and pext- Unfortunately, this may not be possible even in concrete examples. Still, one 
can extract all the information from the functional form of J(X,p ext ). For simplicity let us 
concentrate on the X-dependence and let us consider a generic form of J(X) such as depicted 
in Fig. 1. The attractors so far discussed correspond to the cuts of J with the X axis. 

It is easy to check that the normalization factor for the perturbations, D, introduced in 
Eq. (3.13) in the shift symmetric case is nothing but 



D 



1 



<pK\ dJ(4>, p cx _ t ) 



H 



Pext=const 



(4.14) 



Therefore, we can readily identify non-ghosty regions in the plot of J(4>) as those where 

sgn(D) = a sgn(if) sgn( = +1 , (4.15) 

Therefore, given a fixed choice of branch, a, the perturbations will alternate between ghosty 
and healthy whenever H or the slope J ^ change sign. The roots of J ^ = separate the basins 
of attraction of the different attractors, representing edges between them where D = 0. This 
implies the presence of a pressure singularity which the dynamics cannot penetrate and that 
the sound speed diverges (c^ ~ 1/D for D — > 0). 



The inclusion of nontrivial matter in the discussion is straightforward, it only adds a p C xt direction to the 
plot of Fig. 1. 



- 19 - 




Figure 1: The shift charge J as a function of X for an illustrative shift-symmetric model with 
kinetic braiding; we have suppressed the additional direction p CX f The loci of zeroes of J are the 
attractors: evolution proceeds toward the locus marked as X*. Typically, the attractors alternate 
between being stable and unstable. The basins of attraction are separated by a line of extrema of 
J(X) (marked with the dashed vertical red line). There, the normalization of perturbations vanishes, 
D = Ji(H — 4>k)/H = 0. This leads to a pressure singularity, which the dynamics cannot penetrate, 
and a divergent sound speed, oc 1/D (see the discussion around Eq. (3.16)). The region on the 
left, marked with the jagged line, has the wrong sign of kinetic term, D < and therefore is docs not 
represent a healthy background. 



In the following, we shall concentrate on the evolution close enough to a J = attractor. 
Our approximation requires that two conditions be met: firstly, e J /e* <C 1, where we are 
using the splitting of the fluid introduced in Eq. (4.4); secondly, we must be close enough 
to the attractor such that we are away from any extrema in J, where the backwards-in-time 
evolution would meet a pressure singularity. 

We will therefore study the equation of state of the fluid representing the excess over 
the attractor, with energy density e J = e — e* . Assuming that one can invert J(X,p ext ) 
to find X(J, p ex t), we can simply expand the total energy density e(J, p Gxt ) = (p(J,p cxt )J — 
K(X(J, Pext)) f° r small J to find, at lowest order in J, 

J -I t v, -1 K* x IH-4>*kA 
e J ~ J , where a = 1 - -jy- I — I . 

Notice that the factor in front of J is a function of p cx t, implying that e J does not dilute 
like dust as could have been naively expected and as is seen in the case of ghost condensates. 
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Rather, the equation of state of the off-attractor component is found to be 



w 



J 



(1 + Wext) 



2D* 



K x 



^ext 1 + 2 



dlnE 
dlnX , 



) 



(4.16) 



where we have defined the contribution to total energy density of the external matter as 
^ext — Pext/3-ff 2 . Hence, we see that the off-attractor contribution is also sensitive to the 
external matter, though differently from the attractor. In particular, we see that J behaves 
like dust only when external matter is either absent or when it is just a cosmological constant. 

It is interesting to note that, assuming the H-derivative in the parentheses of Eq. (4.16) 
is not negative enough to change the sign, the energy density stored in the current will 
redshift away more slowly than the external matter. This is irrelevant for the purpose of the 
discussion of the cosmological dynamics, since our approximation assumed that the energy 
density stored in the imperfect scalar will be dominated by its attractor part. However, it is 
true that the minimum of the ratio e J / p e xt will be reached around the time of the transition 
from the J-dominated to attractor-dominated behaviour for the scalar fluid. 

Finally, this discussion has only covered the case when e J is small (as compared to 
either e* or /0 ext ). The opposite regime, when e is dominated by the shift current will have 
significantly different behaviour. The nonlinearity of the equations make a general analysis 
somewhat difficult, so we defer an explicit analysis of this regime to the example of §5. 

4.3 Stability on Attractor 

Specialising to the shift-symmetric case allows us to concentrate on the attractor, where 
significant simplification occurs since J = 0. In particular, using Eq. (4.9), the speed of 
sound can be expressed as 



The second term results from the coupling of the scalar to the other matter content and will 
dilute away as the scalar becomes the dominant source of energy density. When the scalar 
on the attractor dominates the energy density, we can express the sound speed positivity 
condition as 

2 K X *{K*-X*K X *) 

which combined with Eq. (4.12) implies that only one of the branch of the Friedman equation, 
that for which a = sgn</>*, is stable. This equation also displays another difference between 
this theory and ghost condensation: even on the pure-de Sitter condensates the sound speed 
in the braided model is nonzero, and can be adjusted to be positive. The observation that 
an operator like G can lead to > in an expanding universe was also seen in [90] in the 
effective-field-theory language. It should be noted that the linear model presented in section 
§5 does have a vanishing sound speed on the de-Sitter attractor once all the external matter 
dilutes away. In general, this will not be the case, however. 
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5. Example: Imperfect Dark Energy 



In this section, we will describe the properties of the cosmological solutions of the simplest 
model with kinetic braiding which exhibits the interesting behaviour described up to this 
point. We will choose functions K and G to be shift symmetric and linear in X and choose 
the background Minkowski space (X = 0) to be ghosty in order to make the de-Sitter attractor 
the final and stable point of the evolution. The simplest such choice is 



where /i is a coefficient with mass dimension —3. 

Beneath we will in turn discuss the phase space of this theory and the regions of interest, 
the behaviour of the attractor solution, where we will find that wx < — 1 while the sound 
speed squared remains positive. We will then detail how this attractor is approached and the 
behaviour of the energy density of the scalar during this time. Finally, we will discuss the 
regions of parameter space for this model which are not excluded by observations. We invite 
the reader to familiarise themselves with Fig. 2 and its caption: it provides an overview of 
the dynamics of the system which are described in detail in the rest of this section. 

We have named this model Imperfect Dark Energy. The background EMT for cosmologi- 
cal solution has, of course, a perfect-fluid form. However, the EMT for the solution including 
small perturbations can no longer be so described; the fluid is imperfect, as we have shown in 
section §2.2. It is exactly this modification of the fluid's properties which allows this model 
to circumvent certain no-go theorems and stably violate the Null-Energy Condition. 

Other authors have also considered phenomenological models of dark energy which in- 
clude deviations of the energy-momentum tensor away from the perfect-fluid form such as 
anisotropic stress [92, 93, 94]. 

5.1 Phase Space 

The phase-space plot for (ft for an FRW cosmology, extended to show J and H is presented 
in figure 3. We will argue in this section that only the region of (ft to the right of the fixed 
point eft*, where J is positive, is healthy and relevant for the discussion of DE. 
Figure 3 was obtained assuming the branch choice of 



K = -X 




G = \jlX 





J = (ft(3^H(ft - 1) . 



(5.3) 



Choosing the other branch is equivalent to mapping (ft — > —(ft, H — > —H and J — > — J. There- 
fore, for every expanding cosmology on this phase plot, there exists an equivalent contracting 
one with the opposite velocity of (ft. 




2jdX 

Figure 2: The phase space for the Imperfect Dark Energy model Eq. (5.1) in the presence of external 
dust. The red line represents the dynamical attractor (where the shift charge J = 0), with evolution 
occurring toward the final de-Sitter state in the bottom right. The inner purple-coloured region 
surrounding the attractor corresponds to phantom dynamics. Wx < — 1- Outside it, on both sides, 
lie the lighter-coloured non-phantom regions, wx > — 1- The black lines depict example trajectories: 
starting with non-zero shift charge, the system evolves toward the attractor. Should the initial J be 
large enough, the equation of state is non-phantom: the example trajectory on the right crosses the 
phantom divide at the position of the red blob. In this paper, we focus on the region to the right of 
the attractor since the energy density is always positive there. In the region to the left, the pressure 
is negative throughout and it is the fact that the energy density becomes negative in the region e < 
that makes the equation of state non-phantom. Further to the left lies the D < region separated by 
the dashed blue line on which D vanishes and which is the locus of the minimum of J. This blue line is 
a pressure singularity and the region D < is not accessible dynamically from the healthy part of the 
phase space. In the whole region to the right of this dashed blue line, is positive and perturbations 
have positive kinetic term. No solutions exist in the white "dynamically inaccessible" region. 



The evolution of the background will proceed toward one of the two attractors which can 
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be located by solving Eq. (5.3) for J = 0. We obtain for their positions 

ty*Hii = 1 X* = (l^tf 2 )- 1 , (5.4) 
and = 0. (5.5) 

On this attractor the Friedmann equation can be written in terms of the external energy 
density only 

i / r, 2 



H 2 = - ^cxt + yp c 2 xt +3^ 2 J • (5-6) 

An observer who is not aware of the existence of <f> and the attractor would be highly confused 
by this modification of the cosmological dynamics. 

Since J must approach zero in an expanding universe, the basins of attraction will be 
delimited by the extrema of J(4>). The system cannot penetrate these separatrix lines in 
the phase space since pressure is singular there: as already discussed in section §4.2, the 
normalisation of the kinetic term for perturbations, D, is proportional to J^/H and vanishes 
on these extrema, causing both <j> and the sound speed to diverge. Also, this implies that 
the basin of attraction around <j) = 0, since J has a negative gradient there, has ghosty 
perturbations. 

For all the values of <fi < (QfiH)" 1 (to the left of the blue line in Fig. 3), the energy 
density of the imperfect scalar is negative and we will not consider these phase-space regions 
further, since this cannot give the desired dark-energy phenomenology. The remainder of the 
negative- J region does have positive energy density but the energy that can be stored in the 
scalar in the far past can only be very small and to all intents and purposes the dynamics 
would have been indistinguishable from the attractor for all of the observable history of the 
universe. 

Hence, our discussion is going to concentrate on the region (f) ^ 4>* where both J and the 
scalar energy density, e, are positive, the perturbations are not ghostly and the background 
evolution in an expanding universe occurs toward 0* from above. We will show that in this 
region there are no gradient instabilities and the scalar field's evolution provides a viable and 
interesting model for dark energy. 

5.2 Attractor Behaviour 

In this section, we will assume that the initial value of J was small enough such that the 
behaviour of the scalar at times relevant for our observations of the cosmology is effectively 
that of the attractor. As discussed in section §5.1, we will focus solely on the healthy attractor 
at 0*. 

On the attractor, since J = 0, the Friedmann equation simplifies. From Eq. (4.10), 

3H 2 = A* + p cxt , (5.7) 
allowing us to define the cosmological quantities 

= 3^ = ^tf 4 )- 1 n cxt = , (5.8) 
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Figure 3: Plot of the linear model's shift-charge density J (black solid line) and the Hubble parameter 
H (red dashed line) as a function of </>, for constant p e xt- H is negative to the right of the red vertical 
line while the energy density in the scalar is negative in the whole region to the left of the blue vertical 
line. The jagged line marks the region where the perturbations are ghosty. Since D oc J^H^ 1 , there 
are only two ghost-free regions: that around 0* with > and H > and that to the right of the 
diagram with J \ < and H < 0. 



representing the contributions to the total energy density of the imperfect scalar and external 
matter, respectively. 

The Friedmann equation Eq. (5.2) implies that when the Universe empties and p cxt — > 0, 
the requirement to keep the content of the square root positive bounds X from below, 



this puts a lower limit on H, i.e 



X*>±. (5.9) 

fl2> 3^' (5 ' 10) 

The somewhat counterintuitive result is that a small parameter /x, i.e. if the G(X)\I\(f> operator 
is suppressed by a large mass scale, results in a universe with a large effective cosmological 
constant and therefore is excluded by observation. The mass scale suppressing the new term 



must in fact be small 19 . The only way of evading this constraint is by introducing an external 
cosmological constant which will relax the bound, also from below. If there is no external 
cosmological constant, and the current acceleration of the universe is being driven by an 
attractor such as the one discussed then the energy scale suppressing the G\3(j> term is 

/.T 1 / 3 ~ (iZ^Mpi) 1 / 3 ~ l(T 13 eV . (5.11) 

Despite being small this mass scale has the benefit of being technically natural: the radiative 
corrections to this term should be small since it is an irrelevant operator, as opposed to 
quintessence where the scalar mass suffers from quadratic divergences. 

We can explicitly confirm that the perturbations on the attractor have a positive kinetic 
term. Evaluating Eq. (3.13) for the linear model at we obtain 

D ' = 1+ 54^ = 1 + !J *' (5 - 12 > 

Now we turn to the study of the equation of state of the imperfect scalar. The assumption 
that X moves on the attractor implies that 

X* = ((6/i 2 -ff 2 ) -1 ) = -2X*—. (5.13) 

H 

The presence of the coupling to the Ricci tensor in the scalar equation of motion Eq. (2.5) 
resulting from gravity braiding implies that the scalar is sensitive to the presence of external 
matter. This is manifested by the dependence of the position of the attractor on the evolution 
of the Hubble parameter. Differentiating Eq. (5.7) with respect to time leads to 

_ = -- + w ext ) - n xj ^ . (5.14) 

Solving for H/H 2 , we can write down the equation-of-state parameter for the scalar: 

1 + W * x = -(1 + Wext )^^. (5.15) 
1 + U x 

which matches the result we could have obtained from Eq. (4.8). In particular, we find the 
following limits 

1 + w* x » -(1 + w ext ) fi*<l (5.16) 

l + W^-ll + WatlyWO ^cxt«l (5.17) 

As we discussed in section §4.1, while the imperfect scalar is on its attractor, it behaves as 
a phantom and its energy density grows with time. In particular, while subdominant, it 



19 It is of course quite possible to have a very large mass scale suppressing GO</> provided that the final 
vacuum is the trival Minkowski one and dark energy is a result of a cosmological constant 
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anti-tracks the equation of state of the dominant energy source, by which we mean that its 
equation of state is related to that of the external matter by Eq. (5.16). As the dark energy 
begins to dominate the energy budget, it approaches the de-Sitter equation of state from 
below, where it remains as a pure de-Sitter condensate. 

It is important to stress that while this behaviour is phantom, the perturbations to the 
fluid with kinetic braiding are not ghosty, by virtue of Eq. (5.12). In addition, the fluid is free 
of gradient instabilities: the speed of sound squared is positive. Using the results of section 
§4.3, we obtain for the sound speed on the attractor 

with the following limits in the early and late universe 

<&. = |+u>«t nx<l, (5.19) 

cl = f | + u; ext) ^cxt « 1 • (5.20) 

The speed of sound at the attractor is always positive and is determined by the external 
energy density. 

The history of evolution of the dark energy is presented in figure 4 while the hydrody- 
namical properties of the attractor of the fluid of the imperfect scalar are summarised in 
table 1. 



w x 








n x ~ o 


n x ~ i 


MD: 


-2 


— 1 — O ex t/2 


RD: 


-7/3 


1 2 o 


QdS: 


-1 - (1 + W ext ) 


-1 - |(1 + Wext)^cxt 


c 2 


MD: 


4/3 


5 O 

12 ext 


RD: 


5/3 


1 o 

2 11 L cxt 


QdS: 


| + (1 + Wext) 


(2 + 3(1 + w ext )) 



Table 1: Limiting values of the equation of state parameter and the speed of sound for the scalar 
evolving on its attractor. MD and RD refer to the domination of the external fluid by matter or 
radiation, respectively, QdS is the inflationary quasi-dc-Sitter phase. 

From the point of view of structure formation, the above results show how fundamentally 
different the gravity-braided fluid is to a cosmological constant: once the contribution of the 
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Figure 4: Evolution of dark energy properties in a universe containing dust, radiation and the linear 
imperfect scalar evolving on its attractor throughout the presented period. During matter domination 
wx = —2, while wx = —7/3 during radiation domination. The sound speed is superluminal when 
the scalar energy density is subdominant, becoming subluminal when fix ~ 0.1 and wx ~ —1.4. It's 
interesting to note that a model [36] based on higher-order covariant galileon terms exhibits the same 
behaviour of the its equation of state as Imperfect Dark Energy. 



external matter decreases significantly below f2 ext = 1, the sound speed of the fluid will 
begin to decrease. The perturbations begin to propagate subluminally once fix > 0.1. As 
the fluid evolves towards its final de-Sitter state, the speed of sound decreases allowing it to 
collapse and cluster. This clustering of the dark-energy fluid should enhance the growth of 
structure significantly more than is observed in standard quintessence scenarios where the 
sound speed is always c|? = 1 (see, for example Ref. [95] for detailed discussion of the impact 
of non-canonical sound speeds on observations). 

5.3 Approach to Attractor 

It is key to consider precisely the manner in which the attractor described in the previous 
section is approached during the evolution of the cosmology. We will evaluate the hydrody- 
namical parameters for the imperfect scalar whose energy density is dominated by the shift 
current and show that it behaves in such a way so as to rapidly dilute during a period of 
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inflation. However, should a significant shift-current density remain after inflation, we will 
show that it redshifts more slowly than radiation, achieving a local maximum in energy den- 
sity at matter-radiation equality, acting as a model of early dark energy. Then the energy 
density in the scalar dilutes with respect to the dust reappearing once the scalar evolves to 
the vicinity of the attractor and attaches to the solution described in section §5.2. We shall 
refer to the period during which the scalar's dynamics are dominated by the shift-current 
density contribution as J-domination. 

For J / and > 0* , it is useful to define a dimensionless charge parameter 

Q = 6fiHJ. (5.21) 

Our whole region of interest is then spanned by Q ^ 0. It turns out that the transition from 
J-dominated to attractor-like behaviours occurs when Q ~ 1. 

We can re-express all the quantities using Q and in particular rewrite all the dimensionless 
parameters of the cosmology using just Q and Qx- Then, large-Q limits can be found which 
represent the leading-order behaviour of J-dominated scalars. For example, 

^=6^(1 + VT+2Q)-|§, (5.22) 
( l + 2Q + y/l + 2Q(1 + Q)) ~ y^jl , (5.23) 



<Q 



(5.24) 



The above expression for D can be clearly seen to be always positive which verifies explicitly 
that the perturbations are healthy in the whole of the phase space of interest. 

Such manipulation makes it possible to calculate the hydrodynamic parameters of the 
scalar fluid during J-domination . In particular 

1 -w cxt (l -Six) / KOK x 

w x = A + nx ' (5 - 25) 



2 20 + i2 Wext (i - n x ) - (4 - n x )n x 



m + n x ) 2 



(5.26) 



The various limits of the above expressions are presented in table 2. Most importantly, just as 
in the case of the attractor, c 2 > for all the phases of evolution of the cosmology, confirming 
that the perturbations of this fluid do not exhibit gradient instabilities. It can be seen from 
the expression above that the speed of sound is subluminal when the energy in the scalar field 
is J-dominated. 

Another important feature of this model is that fact that it is that the scalar is able to 
cross wx = — 1: a J-dominated fluid has wx > 0; as the attractor is approached, the scalar 
transitions to evolve with a phantom equation of state. Finally it approaches wx = — 1 from 
below. 
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RD: 
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1/4 
1/6 

1 (l- I(l+™ ext )) 




g (1 - gfi ext ) 
5 (l — f^ext) 
" (Wext - g)^ext) 


c 2 


MD: 
RD: 
QdS: 


5/12 
1/2 

± + ±(l + ™ ext ) f- 


r 125 


23 12f2 cx t 
75 125 

23 56Q cxt 
75 375 

(2 + 5(l + «»«t)) 



Table 2: Values of the equation-of-state parameter and the speed of sound for the off-attractor scalar 
dominated by the shift-current J. MD and RD refer to the domination of the external fluid by matter 
or radiation, respectively, QdS is the inflationary quasi-de-Sitter phase. 



We now turn to the discussion of the approach to the attractor in the early universe. We 
assume that the initial conditions were such that the innaton dominated the energy density 
sufficiently to be able to put the universe in a quasi-de-Sitter inflationary phase. However, 
we will also take the initial energy density stored in the imperfect scalar to be similar to that 
of the innaton, since these seem to be the most natural initial conditions. 

Given the model parameters required by the observation of dark energy today, this re- 
quires that most of the energy density initially be stored in the shift- current; the contribution 
of the attractor's energy density during inflation is tiny (since it's inversely proportional to 
H 2 , see Eq. (5.4)). We will now show that a short period of inflation is sufficient to dilute 
this energy density to the extent that the dynamics of the imperfect scalar is effectively that 
of the attractor (i.e. Q < 1) already significantly prior to the acceleration era. 

Using the large-Q expression Eq. (5.23) and dropping numerical coefficients gives 



ff 4 Q 3 / 2 
H 2 



e ~ . (5.27) 



Then, since pi n fi ~ H 2 ni{ , the Hubble parameter during inflation, we have for the initial value 
of Q 

Qi„it~(tt!?) 2/3 (§^) 8/3 (5.28) 

During the quasi-de-Sitter phase of inflation, H remains approximately constant, so Q = 
6fiHJ oc a -3 . Therefore, given J\f e- folds of inflation, we have for Q at reheating 

QRH-Qinite" 3 ^ (5.29) 
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Assuming that reheating occurs without a significant change in H, we can then write down 
the value of Q at matter-radiation equality, given that during radiation domination Q oc 



Taking O^* ~ 1 and H inR ~ 10 16 GeV implies that already for N ~ 17 is Q eq ~ 1. As 
a result, under the usual M ~ 60 e-folds of inflation, the scalar approaches the attractor 
significantly before even BBN. 

It would seem natural that under any such initial conditions as those discussed above, we 
should just take the scalar as being on its attractor throughout the whole observable period in 
the evolution of the universe, and therefore for its energy density to be absolutely negligible 
until very recently. However, since the initial conditions are unknown, we will discuss the 
potential effects of non-zero J in the recent universe. 

After reheating, the scalar redshifts with wx = 1/6, i.e. more slowly than radiation. 
During matter domination this becomes wx = 1 /4. There is therefore a local maximum in Qx 
at matter-radiation equality, provided that the the scalar energy density is still J-dominated. 
This behaviour allows there to be a significant £lx at equality while evading the constraints 
arising from the rate of expansion during BBN. Such a contribution to early dark energy 
("EDE") could potentially be measurable: it would affect the size of the acoustic horizon for 
both the cosmic microwave background and baryon acoustic oscillations, and alter the rate of 
growth of structure at redshifts close to matter-radiation equality and substantially change 
predictions (see for example [96]). A typical evolution of the parameters for the imperfect 
scalar under such conditions is presented in figure 5. 

5.4 Observationally Viable Parameter Space 

The linear model proposed in this section has two parameters: the value of \x and the amount 
of shift charge remaining today. The model can then be reparameterised in terms of the 
dimensionless £lxo arid Qq, the values of those parameters today. 

A large Qq is clearly excluded, since that would give wx > and therefore dark energy 
would dominate the energy density in the past. Realistic cosmologies have dark energy 
dominating the energy density today with small values of Qq, such that the evolution is on 
the attractor in the recent past. 

Beyond the redshift of a few, the evolution is allowed to be off-attractor, provided that 
the energy density Qx at equality is not excessive (Cl e ^ < 0.1, say 20 ). This limits the sort of 
cosmologies that are permitted to be ones described by a narrow range of parameters wxo 
and w a Q. The model discussed here is not a constant w a model. However, for a convenient 
translation of the more frequently used parameter spaces, we include a plot of approximations 
to the above parameters for viable cosmologies in figures 6 and 7. 

20 This is of the order of magnitude of the constraints given in [97, 98, 99, 100], but clearly the actual 
constraints would need to be recalculated in a full analysis. 



a' 5 oc H 5 / 2 




(5.30) 
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Figure 5: Evolution of dark energy properties in a universe containing dust, radiation and the linear 
imperfect scalar. The energy density in the scalar is J-dominated until a transition during the matter 
domination epoch. This allows the scalar to increase its contribution to the total energy budget 
throughout radiation domination (wx = 1/6) and provide an early dark energy peaked at matter- 
radiation equality, from whence it begins to decline with wx = 1/4. The transition to the attractor 
behaviour is rapid (the time at which the transition occurs depends on the initial value of J). The 
equation of state crosses wx — — 1 and the scalar energy density begins to grow. The final stages of 
evolution are on the attractor and are effectively identical to those presented in figure 4. The speed 
of sound remains positive throughout. 

Note added: Shortly after the publication of the pre-print of this paper on the arXiv, 
Ref. [102] appeared which discusses the model (1.1) in the context of inflation. In this ref- 
erence, the authors present the quadratic action for perturbations and find the conditions 
under which the perturbations remain healthy. These results agree with ours. 
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Figure 6: Plot of range equation-of-state parameters allowed by requiring that the cosmological 
model be viable: Models are chosen by requiring that 0.1 < Q m < 0.5 and Slxcq < 0.1. The shading 
contours correspond to the energy density of dark energy today flxo- Two parameterisations of dark 
energy behaviour are shown: wx and w' x evaluated today, which is a reasonable approximation in the 
region of redshift to which supernovae are sensitive, and wx evaluated today and at z — 1/2, which is 
an alternative parameterisaton [101]. The requirement that the energy density in the imperfect scalar 
at matter-radiation equality be small, ft c x < 0.1 forces the value of the shift charge to be small today 
Qo < 10~ 2 . This means that in the most recent history, the evolution has effectively been on attractor 
or very close to it and the permitted value of wx is very restricted and determined to all intents and 
purposes by fl x . 
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A. Action for Cosmological Perturbations 

In this section, we derive the quadratic action for perturbations of the full scalar-tensor 
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Figure 7: Plot as in figure 6, but with shading representing the contribution of the imperfect scalar 
to energy density at matter-radiation equality. We choose to cut the parameters such that the con- 
tribution to this early dark energy at that time is no larger than 10%. It can clearly be seen that 
values of wx closer to —1 arc obtained when the shift charge is larger, but this leads to more early 
dark energy, eventually disagreeing with current constraints. 



system around an arbitrary homogeneous and isotropic background. We follow the method 
introduced in [103]. 

Let us start by writing the action Eq. (2.1) in a slightly more convenient form. Integrating 
by parts twice the G B term, one obtains 

S=\f d A x^ i ] R + 2K + 26GV2X + G ^ } , (A.1) 

where as before 6 = V^u^ is the trace of the extrinsic curvature of the slicing (j) = const with 
the normal = V ^(p/ y/2X . 

The computation of the quadratic action is most convenient in ADM variables, where 
the metric is written as 

As 2 = N 2 dt 2 - hijirfdt + dx l )(N j dt + dx j ) (A.2) 
where is the lapse function, iV* is the shift vector and hij is the induced metric on the 
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slices of constant time t. The extrinsic curvature of these slices is Eij/N with 

= -hij - D { iNj) 

where Di is the covariant derivative associated with hij and the indices are raised and lowered 
with h^. 

The usefulness of these variables is that in the 'unitary' gauge, defined as <p = (f>o(t) (that 
is 8(f) = 0), the slicing defined by t and <p coincide. Consequently, the i,j components of the 
spatial projector and of the extrinsic curvature K^ v =_L? V p u a coincide with hij 
and Eij/N respectively. Furthermore, in this gauge, the full unperturbed form of X is simply 

X = 4>1/{2N 2 ) , 

which depends only on 5N. 

The background solution is given by hij = a 2 (t)5ij, N = 1, N % = and cp = ^o(*)- It is 
also convenient to decompose the spatial metric as 

h ■ ■ — p 2a(t)+2((t,x) £ 

where a = loga(i) and det(/ijj) = 1. The variable C is the curvature perturbation. 
To fix the remaining gauge freedom it is useful to impose 

d% 3 = , 

so that hij is a 'transverse-traceless' variable, describing the tensor modes. 
The action in ADM variables is 

S= l -f d 3 xdtv^ |iV(i? 3 + 2K(X)) + 1 [El - E 2 ) + (A.3) 

+2GV2XE + ^{X4> + N l diX^j I , 

where R3 is the Ricci scalar for hij, we have dropped terms containing di4> which is always 
vanishing in the chosen gauge and we remind that X = X(N). 

The strategy devised in Ref. [103] to obtain the quadratic action was to first compute 
8N and 5N l = N l at first order in terms of £ and hij using the Hamiltonian and momentum 
constraints and then to substitute them in the action and expand it to quadratic order. The 
expression of 5N and N l to second order is not necessary because those terms are multiplied 
by the constraints at zeroth order, which identically vanish. 

The momentum constraint is 

Di ^-^ +GV ^ hii)j __G_ gtX . =0 

It is convenient to separate the shift vector into its transverse and longitudinal parts, iVj = 
Nj -\-diip with D l Nf = 0. Expanding the momentum constraint at linear order, its transverse 
and longitudinal parts give 

Nf = (A.4) 
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and 



(H-G X XV2X)\ SN = (. 



(A.5) 



respectively, tp does not appear in the quadratic action since it enters the action (A. 3) 
through N l , which only appears (up to total derivatives) in the combination N' 1 times the 
linearised momentum constraint, which is set to vanish. As usual, the subscript o denotes 
the background value, e.g., Xq = 4>\j2 etc. Notice here that the relationship between the 
perturbations of the lapse and the curvature, Eq. (A.5), is changed from the one in It-essence. 
One could expect that this will impact the normalisations of inflationary perturbations in a 
new way. 

Substituting the momentum constraint into the action, dropping terms which will only 
contribute at orders above second and integrating by parts we obtain the partially reduced 
action 

S=^JVh N(R 3 + 2K) + ^ (A% - A 2 ) - 2&V2X . (A.6) 
where Ay = hij/2. Substituting for perturbed variables we find: 

t e 3(ff+c) { (1 + SN) (2K(<f>, X(l + 5N))e~ 2( - a+0 

(hf 



S 



4AC + 2(dC) 2 + 



{dhf 



+ 



1 

1 + 6N 



6{H + C) 2 + - 2G((f>, X(l + 5N))^/2X(1 + 5N)} . (A.7) 



Expanding to second order, after some manipulations and using the background equations 
of motion, one finds the quadratic action 



5(2) 



c 2 



a 



+ 



1 



: (d k h 



I.) I 



where the normalization of the kinetic term is 
2XD 

A 



{H-^XGx? 

D = K X + 2XK XX - 2G^ - 2XG x<t> + ^H(G X + XG XX ) + 6X 2 G 2 X 

and the sound speed for the scalar perturbations is 

2 _ <j>XG x (H - <j)XG x ) -{H- <])XG x y 
° s XD 



(A.8) 

(A.9) 
(A.10) 

(A.ll) 



where it is understood that everything is evaluated on the background solution. The Fried- 
mann equations then have to be used to eliminate H, remembering that this derivation was 
performed under the assumption of the imperfect scalar's providing the total energy content 
of the cosmology. In this way the above reduces to 



K x - 2G + 2XG^ X + 24>{G x + XG XX ) + 4<t>HG x - 2X 2 G 2 X 



D 



(A.12) 
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The dependence of the sound speed on cj) remains and the equation of motion for the imperfect 
scalar has to be used in order to complete the evaluation. Generically, therefore, the sound 
speed will depend on the dynamics of the whole system and is not just a property of the fluid. 

Finally, note that the properties of the tensor modes are not affected. The sound speed 
is still 1, as expected because we are only modifying gravity through the trace of the extrinsic 
curvature, which is independent of hij. 
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